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In this work a supersymmetric cosmological model is analyzed in which we consider a general 
superfield action of a homogeneous scalar field supermultiplet interacting with the scale factor in a 
supersymmetric FRW model. There appear fermionic superpartners associated with both the scale 
factor and the scalar field, and classical equations of motion are obtained from the super-Wheeler- 
DeWitt equation through the usual WKB method. The resulting supersymmetric Einstein-Klein- 
Gordon equations contain extra radiation and stiff matter terms, and we study their solutions in 
flat space for different scalar field potentials. The solutions are compared to the standard case, in 
particular those corresponding to the exponential potential, and their implications for the dynamics 
of the early Universe are discussed in turn. 
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I. INTRODUCTION 

From the latest observations ,we do know that about 
95% of matter in the Universe [1| is of non-baryonic na- 
ture, and the rest is constituted by radiation, baryons, 
neutrinos, and all other particles we understand well in 
the Standard Model of Particle Physics. The most suc- 
cessful models until now is the so called Lambda Cold 
Dark Matter (ACDM) modeljj^, which is able to ex- 
plain and to fit reasonably well all cosmological observa- 
tions. 

In the last decades cosmologists have made use of 
scalar fields in the description of various aspects of cos- 
mology. They are known in models of inflation, and more 
recently in models of dark energy, see[8, 9] and references 
therein. But scalar fields have been considered too for 
models of dark matter [Tot [Tl|. 

This flexibility of scalar field models to describe dif- 
ferent phenomena comes from the properties of the self- 
interacting scalar field potential V((p) that is specified 
for each model. Currently, there is no underlying prin- 
ciple that uniquely specifies the potential for the scalar 
field and many proposals have been considered [cl H-1 [l3j . 
Some were based in new particle physics and gravita- 
tional theories, other were postulated ad-hoc to obtain 
the desired evolution. 

On the other hand, the physics required to understand 
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the early Universe should be necessarily rooted in a the- 
ory of quantum gravity. Futhermore, it would probably 
be adequate to consider scenarios where both bosonic 
and fermionic matter fields would be present on an equal 
footing. 

In considering the quantum creation of the Universe 
we are of course dealing with the earliest epochs of the 
Universe's existence, at which time it is believed that 
supersymmetry would not yet be broken. The inclusion 
of supersymmetry could therefore be vital from the point 
of view of physical consistency. 

For these and other physical reasons supersymmetric 
quantum cosmology emerged as an active area of re- 
search. The first model proposed [___| was based on the 
fact that, shortly after the invention of supergravity, it 
was shown [ill flo| that this theory provides a natural 
classical square root equations and their corresponding 
Hamiltonians. 

A second method later proposed was a superfield for- 
mulation, in which is possible to obtain the correspond- 
ing fermionic partners and also being able to incorpo- 
rate matter in a simpler way[l__-[___|. A third method 
allows us to define a square root of the potential, in 
the minisuperspace, of the cosmological model of inter- 
est and consequently operators which square results is 
the Hamiltonian [20| - K22| . 

So, in the same way that we seek a desirable scalar 
field potential to explain the evolution (and early times) 
of the Universe from the point of view of standard Gen- 
eral Relativity, we can reconcile these requirements along 
with the ideas of local supersymmetry using now super- 
potentials. For this purpose we need to model a super- 
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symmetric quantum cosmological landscape and see what 
happens now with the super expansion factor and with 
the super-scalar fields. It is then important to find the 
influence of the " fermionic" variables in these supcrficlds 
and how they would alter the dynamics of usual cosmo- 
logical models. 

In this work, we consider a Hamiltonian for a ho- 
mogeneous super scale factor, which is a supermultiplet 
(with four components and different signs) in supergrav- 
ity N = 2, and that interacts with a super scalar field 
(also a supermultiplet) [19J . We shall promote this Hamil- 
tonian to be an operator, representing the Grassmann 
variables by matrices, then by means of the WKB proce- 
dure we find two (independent) classical evolution equa- 
tions. Those associated with the scalar field are obtained 
through Hamilton equations. 

This procedure gives us a modified Einstein-Klein- 
Gordon (EKG) set of equations (that we call SUSY-EKG 
equations) due to the indirect presence of the "graviti- 
nos" and the "fermionic" variables corresponding to the 
scalar field, which are inherently contained in each entry 
of the supermultiplets. From a phenomcnological point 
of view, the new extra terms in the model offer different 
kind of components that behave as radiation and stiff 
matter. 

The paper is organized as follows. In Sec. lU we out- 
line the procedure that allows us to define the superfields 
associated with the expansion factor and the scalar field 
so that we can generalize the usual action of Cosmology. 
We find the Hamiltonian of the system, which already 
contains extra terms depending upon the Grassmanian 
variables associated with the scale factor and the scalar 
field. 

The Grassmanian variables are represented as matri- 
ces, and then the Hamiltonian operator is a matrix itself 
with four components. We focus our attention in its two 
independent components, and apply to them the usual 
WKB method to get classical equations of motion. Some 
solutions are found for the cases of a free scalar field and 
of a constant scalar field potential which is negative def- 
inite. 

Sec. Mil is dedicated to the analysis of the case in which 
the scalar field is endowed with an exponential super po- 
tential. We first show that there is an exact scaling solu- 
tion, in which all energy terms behave like stiff matter. 
To have a complete picture of the solutions, the equa- 
tions of motion are written as a dynamical system, and 
we study its critical points and general trajectories in the 
phase space of the resulting variables. 

Finally, Sec. IIVI is devoted to conclusions and com- 
ments about the general properties of the SUSY-EKG 
equations and their solutions. 



II. THE SUSY-EKG EQUATIONS FOR A FRW 
UNIVERSE 



In this section, we describe the main features of SGR 
cosmology, and for this we will write the supersymmetric 
version of the WDW equation according to the superfield 
method outlined in the introduction. 



A. Mathematical background 

For a homogeneous and isotropic universe, we write the 
Friedmann-Robertson- Walker metric as (in units with 
c=l), 



ds 2 = -N(t)dt 2 + a 2 (t) 



dr 2 



1 — kr 2 



r 2 dil 2 



(1) 



where a(t) is the (time-dependent) scale factor, N(t) is 
the lapse function, and k is the curvature constant. Then, 
we can write the total action representing a (real) scalar 
field 4> endowed with a scalar field potential V(4>), and 
interacting with the expansion factor as 



S = 



G 



8vrG 



~ + \kNa) l ll s ,(<T». (2) 



The equations of motion arising from this action, for iV 
1, are 



H = -—<?, 



together with the (constraint) Friedmann equation 



H 2 = 



y 2 + V{4>) 



(3a) 
(3b) 

(4) 



where H = a/a is the Hubble parameter, and k 2 = 8irG. 
Eqs. (I5al) . (|3bl) . and are the representative equations 
of motion of a FRW universe driven by a scalar field. 

We want now to review the procedure that one of us 
and collaborators have followed to construct a superfield 
action for the FRW model interacting with a (homoge- 
neous) scalar supermultiplet 19], and from this the su- 
pcrhamiltonian associated with it. 

The most general superfield action [13, EH, HI! has the 
form 



s 



G 



1 , Vk 



A 2 



2 JV 



V^V^ -2A 3 g(<f>) 



dr] dr] dt 
dr] dfj dt , (5) 



where k = 0, 1 denotes flat and closed space, and k 2 = 
8ttGn, where Gn is Newton's gravitational constant. 
The units for the constants and fields in this work are 
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the following: [n 2 ] = £ 2 , [Af] = £°, [A] = I 1 , [$] = £~\ 
[g(&)] — £~ 3 , where £ corresponds to units of length. 
Besides, T> v =8^ + ifjdt and 2?^ = —d^ — irjdt are the su- 
percovariant derivatives of the conformal supersymmetry 
N = 2, which has dimension [D n ] = [V n ] = f 1 ' 2 . 

For the one-dimensional gravity superficld Af(t, ij, ff) 
(Af = Af' ) , we have the following series expansion, 



Af(t, r), ff) = N(t) + ir)ip'(t) + ifjtp'(t) + T]fjV' . 



(6) 



where N(t) is the lapse function, and we have also intro- 
duced the reparametrization ip'(t) = iV 1 / 2 (i)?/'(t), and 
V'(t) = Af(t)V(t) + ij)(t)ip(t). The Taylor series expan- 
sion of the superfield A has a similar form, 

A(t, 7], ff) = a(t) + it]X'(t) + ifj\'(t) + rjfjB' , (7) 

where a(t) is the scale factor, \'(t) = KN 1 / 2 (t)X(t), and 
B'(t) = KN(t)B(t) + (l/2)«(V5(i)A(t) - ip{t)X(t)). Like- 
wise, the scalar superfield rj,fj) may be written as 

$(f, r,, fj) = <j)(t) + irjx'it) + ifj X '(t) + m?' , (8) 
where _x'(i) = AT 1/2 (i)x(t), and P(t) = N{t)F(t) + 

(i/2)(f(t)x(t)-m X (t)). 

As it was shown in Ref.[17J, we now expand the ac- 
tion §5§ in terms of the superfield components ©, Q), 
and ijSj), and integrate over the Grassmann complex co- 
ordinates r\ and fj. If we redefine 

A(t)-> V 1/2 WA(t), xW^«" 3/2 WxW, (9) 

it is possible to find the Lagrangian, and from it the 
superHamiltonian can be constructed, namely, 

nj i/2tt i/2-i-r 3fca 1 r - 1 n 2 

n = --a' n a «/ n a - — --— [a,a] + ^ 



n v ([A, x ] + [A,x]) 



16cr 



[A, A] [x,x] 
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"la 3 

~ x] + yS [A, A] + 67fc 5 fo>) a 2 

. a 3^M + ^ 2 5 (^) [x,x] + ^^[x,x] 



([A,x] -[A,x]) • 



(10) 



2 <9v? 

where the scalar field potential reads 

%)^(|>) 2 - 3K y M . (id 

Notice that, in general, the scalar potential (fTT|) is not 
positive scmi-dcfinitc. The relevant term in Eq. (11 1 1) is 
which is related to the superpotential and whose 
form shall be chosen appropriately for the cosmological 
model under study. 



In the quantum (canonical) formalism the Grassma- 
nian variables A, A, \, and x, — by the anticommutators 
as 



{A,A}=--, {x,x} = l, 



(12) 



and they can be considered as generators of the Clifford 
algebra, as well as the commutators 



[a, n o ] = -i , 



,n ri 



(13) 



We can choose a matrix representation for the 
"fermionic" operators A, A, x, and x, m t ne form of a 
tensorial products of 2 x 2 matrices, 

A=y|o-_®l, A = -w|ct+®1, (14a) 



X = cr.3 ® cr_ , x = ^3 ® 0"4 



(14b) 



where <t± = (01 ± 1172) /2, oi, 02, and 173 are Pauli ma- 
trices. 



B. The classical landscape 

As we have already mentioned in the introduction, our 
objective is now to construct the classical equations that 
corresponds to the Hamiltonian (|10p . 

First, we promote it to an operator W by realizing 
the "fermionic" variables as the matrices (|14l) and, as 
usual, II Q = id a and 11^ = id,/,. By these means, we will 
get a quantum Hamiltonian operator that should fulfill 
T-L\^) = 0. Because the matrices in (|T4l are 4x4, the 
wave function will have four components. 

It can be shown [l9j that the components ^1 and ^4 
satisfy independent equations, whereas ^2 and ^3 ap- 
pear coupled in the other two differential equations. In 
this work we focus our attention in the former case, and 
apply separately to and 'J 4 the WKB method, so that 



(15) 



With this we shall find the classical equations of mo- 
tion associated to the components ^1 and ^4, which also 
correspond to the classical Hamiltonian; from this, classi- 
cal equations can be obtained for the scalar field ip. Thus, 
the classical equations of motion are 

.. dV ,3K 2 dg(v) r Vkdg(if) 
p = —inp — — ± — — — 6- 



dip 2 a 3 dp a dp 

1 d 3 g(p) 



a 3 dp 3 



(16a) 



9 K 2 ip 2 k 2 , N k n 2 \fk ft 4 



6 3 

A 



3a 4 32a 6 



,- 9M+2 ^ 9M± ^^), 16b) 

The upper (lower) sign in Eqs. (|16[) corresponds to the 
(quantum) equation for ^1 (^4). 
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We have now SUSY-EKG classical equations that, due 
to the presence of variables A and A associated with the 
"gravitino" and the "fermionic" variables \ an d X asso- 
ciated with the scalar field, considerably differ from the 
standard EKG equations <j3j> of General Relativity. 

There are extra terms, due to supersymmetry, behav- 
ing like radiation (a~ 4 ) and stiff matter (a -6 ), which 
should be dominant at very early times, whereas other 
terms show a combination of scale factor powers medi- 
ated by the presence of the superpotential g(4>) and its 
derivatives. 



C. Simple classical examples 

As a first instance of a solution, we will consider the 
case g(<p) = 0, which also corresponds to a null potential, 
V(ip) = 0. We also set the curvature term k = 0. The 
solutions of Eqs. (|16|) are 



= <A)(ao/a) 3 , 



a 3 (t) = a 3 , + 3 



2 • 2 fi 4 \ 1/2 

K Z <Pq(Iq X 



G 



32 



(17) 

(t - t ) , (18) 



where to, tpo, and ag are integration constants. The whole 
solution corresponds to stiff matter, and is practically 
the same as in the standard case because the g-terms 
dissappear from everywhere. 

Another more interesting case is that with a constant 
superpotential, g(ip) = go, that corresponds to a constant 
and negative definite scalar field potential, V = — 3ft 2 g 2 . 
The scalar field potential is then an effective cosmological 
constant which is negative definite. 

As in the previous case of the free scalar field, the case 
is simplified because the derivatives of the superpotential 
disappear from the SUSY-KG equation (|16ap . though not 
from the SUSY-Friedmann equation (|16b[) . The cosmo- 
logical solutions can be expressed as 



(p(t) = ipo{a /a) 3 , 



(19a) 



a (*) = r 1 



' 8<P 2 af 
3k 2 



- sin [3K 2 go(t-t )] ^Mb) 



where again ipo and ao are integration constants. 

Because the scale factor is a positive quantity, the only 
acceptable solution is when the amplitude of the sinus 
function is less or equal to one. It is then clear that 
the scale factor has a periodic solution in which ao is 
the amplitude at maximum expansion. In other words, 
Eq. (|19b|) represents an oscillatory Universe. 



III. SUSY COSMOLOGY WITH AN 
EXPONENTIAL SUPERPOTENTIAL 

The possible cosmological roles of exponential po- 
tentials in scalar field models have been thoroughly 



investigated in the specialized literature [2J-|36(, see 
also Hll 13) almost always as a means for driving a 
period of cosmological inflation, but also as possible can- 
didates for dark matter and dark energy. 

Scalar field cosmologies with an exponential potentials 
are, as compared to others, mathematically simple, and 
their solutions have many interesting features. For the 
purposes of this work, we only mention the possibility of 
having inflationary solutions and the appearance of the 
so-called scaling solutions, which are nicely illustrated 
in, for instance, the dynamical system study presented 
in[32j]. 

The inflationary solution for exponential potential is 
the simple power law inflation [l3j], which however never 
ends and needs modifications to provide a graceful exit 
towards a Hot Big Bang model. On the other hand, 
the scaling solution arises whenever the scalar field is 
accompanied by another matter fields, so that both fields 
evolve with a fixed ratio of their energy densities, see for 
instance[H H3, HM1 H . 

In this section we explore in detail the type of solutions 
permitted by our (classical) SUSY cosmological model 
when the scalar field is endowed with an exponential po- 
tential. Our main interest will be to find inflationary 
and scaling solutions. Even though we are not consid- 
ering extra matter fields apart from the scalar field, the 
new terms in Eqs. (|16[) will play the role of companion 
fields which should impose a non-trivial behaviour upon 
the field <p. 

Let us consider the following superpotential and po- 
tential, respectively, 



g( V )=goe- x ^l\ 



V(<p)=V e- x "*, U ^^(A 2 



6 



(20a) 
(20b) 



where the potential parameters were chosen to ease their 
comparison with the standard case; notice that in order 
to avoid a negative definite potential we should impose 
the condition A > v6- The equations of motion (IT6|) with 
an exponential superpotential explicitly read 



a , . , , . \n 3 q 

(p = -3-(p + XkV ± (A - 6) — f 



H 2 



-V 



32a ( 



±(A 2 -6) 



4 

« 9 
12a 3 



(21a) 
(21b) 



A. Exact SUSY scaling solution 

We present here a first (exact) solution of Eqs. (f2"Tj) 
that we shall call scaling solution, because of its resem- 
blance with the scaling behavior exponential potentials 
show in standard cosmology [29j, |33|, |3J|. 

It can be noticed that there is a stiff matter term in 
Eq. (|21bp , and that the superpotential g appears accom- 
panied by factor a~ 3 . Thus, one can foresee that there 
must be a stiff matter solution of the equations of mo- 
tion, so that a ~ t 1 / 3 , as long as g ~ a~ 3 and V ~ a~ 6 . 
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It can be shown, just by direct substitution in Eqs. (l2"TT) . 
that the exact scaling solution is 



where 



a(t) = a (t/t 



,1/3 



mp{t) = - \n(t/t Q ) , 

g(t) = go(t /t) , go = T% 



(22a) 
(22b) 
(22c) 



where ao is an appropriate constant to accomplish 
Eqs. flUDHp. 

The scaling solution corresponds to stiff fluid matter, 
as revealed by the power law behaviour of the scale factor 
in Eq. (|22a[) ; this is probably not surprising, because we 
have already noticed the presence of a stiff-term in the 
SUSY Friedmann equation (|16b[) . This solution is not 
inflationary, but its existence indicates its possible im- 
portance in the early dynamics of the models considered 
here. 



B. Dynamical System Structure 

Our next step is to study the evolution of our SUSY 
classical model in which the scalar field if is endowed with 
an exponential potential. As in the standard case, it is 
possible to perform a dynamical study of the cosmological 
model so that its physically relevant solutions are easily 
unveiled. 

In order to construct a dynamical system for our cos- 
mological model, we follow Ref . [29| . see also[3l|, [13, IH, 
[38j . One first step is to introduce a set of conveniently 
chosen variables which may allow the rewriting of the 
evolution equations as an autonomous phase system sub- 
ject to a constraint arising from the Friedmann equation. 

We choose the following variables, 



Ky3 



which render the Friedmann equation as 



(23) 



F(x, y, z) := x 2 + y 2 + z 2 ± 2^{X 2 - 6)/3yz = 1 . (24) 



The constraint equation ([24| follows from Eq. (|21b| , and 
we see that variable z plays the role of an extra fluid term 
which, contrary to the standard case, see Ref. .29], is not 
trivially coupled to the scalar field variables. 

We shall restrict ourselves to the part of the phase 
space that makes physical sense, and this is the range 
< | a; |, y, z < oo, which is the part that corresponds to 
expanding universes only. Combining expressions (|21[) 
and (|23p , the equations of motion read 



i>x H2 x 



3^ 2 , aVFI ,„ . 
2 — 71 — yZl ^ 



(25b) 
(25c) 



H 
W 



-3x 2 



3z 2 t v^A 2 - 6) 



yz . 



(26) 



Here primes denote derivative with respect to the loga- 
rithm of the scale factor, N = ln(a). The evolution of 
phase space variables x, y, and z takes place only on the 
constraint surface described by Eq. (|2"4"1) . 

Notice that there is a symmetry in Eqs. (1251) with re- 
spect to the double sign (±) and variable y. The case 
with sign (— ) can be obtained from the case with (+) 
if we change y — > —y; and vice versa. As we shall see 
below, we will only study the case with the lower signs 
since only for them is that we can obtain positive results 
for y. 

We show in Tablc|T]thc critical points of the dynamical 
system (|25]) and their (linear) stability properties. 

There are five critical points, in close similarity to the 
standard case, whose main features are described next. 

• Stiff matter domination. The potential variable 
is null, y — 0, and then the dynamical system is 
equivalent to the standard case of stiff fluid mat- 
ter (a~ 6 ) plus a free scalar field ((f) ~ a~ 3 ), so that 
x 2 + z 2 = 1. Particular (well known) cases are: 

— Point A, stiff fluid domination. The scalar 
field variables x and y are both null, and then 
this point represents the complete domination 
of the stiff matter term, z = 1. This solution 
exists also in the standard cosmological case. 

— Points B, kinetic domination. They represent 
the domination of the scalar field's kinetic en- 
ergy, x = ±1, and y = = z. This solution 
exists also in the standard cosmological case. 

— Points C, joint kinetic and stiff domination. 
This is a (restricted) scaling solution of stiff 
nature which is satisfied by all points inside a 
unitary circle on the plane y = 0. This set of 
solutions exists also in the standard cosmolog- 
ical case, but it was missed in the analysis of 
Ref. [291] : it is the line segment x = [—1,1] at 
y = 0. Notice that points A (B) can be seen 
as extreme C-points as z — > 1 (z — > 0). 

• Point D, scalar field domination. It is the coexis- 
tence of the (scalar) kinetic and potential energies, 
x 2 +y 2 = 1, and then the point is located in the uni- 
tary circumpherence on the plane z = 0. Notice, 
however, that the existence of this point requires 
A 2 < 6, which is in contradiction with our earlier 
assumption that A 2 > 6, see the scalar field poten- 
tial defined through Eq. (f20bl) . 

• Point E, scaling solution. This point corresponds to 
the scaling solution in Sec. lIII Al and represents the 
coexistence of all energy terms in the equations of 
motion. It should be noticed that, contrary to the 
present work, in the standard cosmological case the 





FIG. 1. Different views of the phase space of the dynamical system (|25[) for the particular value A = 2.85, so that 6 < A 2 < 9; 
the 3-d surface represents the constraint F(x,y,z) — 1, see Eq. I|24p . whereas the curves are solutions of the dynamical system 
for diverse initial conditions. The dots denote the critical points A, B, C and E shown in Table [I] and the trajectories reveal 
their stability properties as described in the text. In particular, notice the largest (red) loop that encircles point E (see figure 
on the right): it is a homoclinic trajectory that departs from and arrives to the same critical point. 





FIG. 2. Different views of the phase space of the dynamical system (|25[) for the particular value A = 3.5, so that now 9 < A 2 ; 
the 3-d surface represents the constraint F(x, y, z) = 1, see Eq. 1|24[) . whereas the curves are solutions of the dynamical system 
for diverse initial conditions. The dots denote the critical points A, B, and C shown in Table U and the trajectories reveal 
their stability properties as described in the text. Point E does not exist in this case. As in Fig. [I] there also exists the (red) 
homoclinic trajectory that departs from and arrives to the same critical point located at (v6/A, 0, ^/l — 6/A 2 ). However, we 
were not able to show the complete closed trajectory due to numerical limitations. 



scaling solution in the presence of stiff fluid matter 
necessarily requires y = 0, see Table I in Ref.[29j. 



x = x + u, where x = (x, y, z) and u = (Sx, 6y, Sz). The 
equations of motion (|2"5j) can be written as x' = f(x), 



The stability of the points is investigated through lin- 
ear perturbations around the critical values of the form 
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TABLE I. Critical points (a;*, y*, z*) of the dynamical system (|25f) that represents the classical supercosmology of a scalar field 
endowed with an exponential superpotential. See also Fig. [T] for a graphical representation of the phase space. 



Label 


X* 


y* 


2* 


Existence 


Stability 


A 








1 


VA 


Unstable 


B 


-1 








VA 


Unstable 




1 








VA 


Unstable for A 2 < 6 












Saddle for A 2 > 6 


C 







z* 


VA , z,<l 


Unstable 



D 

E 



A/A 



z* 



(A 2 



V1-AV6 
- 6)/[Av^= 



VA, 



< 1 



A^] ^3(A2 - 6)/[AV9^P] 



A 2 < 6 
6 < A 2 < 9 



Unstable for sj\ - 6/A 2 < z < 1 
Saddle for < z < ^\ - 6/A 2 
Stable node 
Stable centre 



which upon linearization reads 



(27) 



where M. is called the linearization matrix. The eigenval- 
ues oj of A4 determine the stability of the critical points, 
whereas the eigenvectors r] of Ai determine the principal 
directions of the perturbations. In general, if Re(w) < 
(Re(u) > 0) the critical point is called stable (unstable). 

In principle, we should study the perturbations of the 
three dynamical variables (x, y, z), but we should remem- 
ber that they are not all independent because they are 
bond together by the Friedmann constraint (|24[) , and the 
same happens for their perturbations. 

The Friedmann constraint defines a two dimensional 
surface upon which lie all physically relevant phase space 
trajectories, and then we will be interested on perturba- 
tions lying also on the constraint surface. In other words, 
perturbations which are perpendicular to the constraint 
surface should be taken away from the analysis. 

We can identify the excluded perturbations by com- 
paring their associated eigenvectors with the gradient of 
the constraint surface at each critical point, 



VF| X =x*i+ (y»± V(A 2 -6)/3z* j 



+ (z* ± ^(A 2 - 6)7%*) 



(28) 



We will only take into account eigenvalues associated to 
eigenvectors for which rj- VF| X = 0. The stability results 
are also summarized in Table [1 

Here we list the eigenvalues of the stability matrix A4 
for each of the critical points, only for the perturbations 
that are compatible with the Friedmann constraint. 

• Point A. It is an unstable point with eigenvalues 



^2 = . 



(29) 



The instability happens only along the eigenvector 
corresponding to cui , which points in the positive y- 
direction. The second eigenvalue is null, and then 
the system is indifferent under perturbations along 
the unitary circumference x 2 + z 2 = 1 . 



• Point B. The stability eigenvalues for the cases a;* 
=Fl are 



Vf(A±V6), 



uj 2 = . 



(30) 



Similarly to the case of point A, only the eigenvec- 
tor corresponding to u)\, which points in the posi- 
tive y-direction, gives information about the stabil- 
ity of the critical points. We notice that the point 
at x = —1 is unstable, whereas that at x = 1 can 
be unstable or saddle, depending upon the value of 
A. Again, the second eigenvalue is null, and then 
the system is indifferent under perturbations along 
the unitary circumference x 2 + z 2 = 1. 

• Points C. The eigenvalues corresponding to = 

tV 1 - z * are 



u>i = 3 ± 



z 2 

- 1 * ) 



^2=0. 



(31) 



These points have stability properties very similar 
to points A and B, but stability also depends upon 
their exact location on the unitary circumference 



special point, corre- 
6/ A 2 and located at 



x 2 + z 2 = 1. 

Notice that there is a 
sponding to Zh = v/l- 
(V6/A, 0, y/l — 6/A 2 ), which marks the instability- 
stability transition of the chain of points C. Because 
of this, there is a particular trajectory that departs 
from and also arrives to the point. This is called 
an homoclinic trajectory, and it is the largest loop 
that encloses point E, see Figs. [1] and [2j 

• Point D. Their stability eigenvalues are 

^i = -i(6-A 2 ), u 2 = - l -{Q-\ 2 ). (32) 

It is a stable point, whenever it exists, i.e. if the 
case A < v6 is allowed. 
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Point E. Its stability eigenvalues are 



6) 



3(A? 

2 V9~~\ 



3/2 



UJ 2 



3(A^ 

2 V9~V 



3/2 



The eigenvalues are purely imaginary for the range 
of existence of the critical point, then it is a stable 
centre. This is confirmed by the closed trajectories 
around the critical point in Fig. [TJ 

The overall conclusion is that only critical points A, B, 
C and E may coexist together in the phase space, because 
the (inflationary) point D is excluded by the (positivity) 
restriction A 2 > 6. What we observe in Figs. Q] and [2] is 
that the ultimate fate for trajectories is to move around 
point E in closed loops, for the case 6 < A 2 < 9, or 
to reach any of the stable points C, for any A 2 > 6. All 
possible solutions in the phase space represent, in general, 
stiff matter solutions. 



IV. CONCLUSIONS 

In this work we have considered a supersymmetric ex- 
tension of the action of general relativity for a scalar field 
interacting with the scale factor of the Universe. For this 
purpose, we have introduced a superfield formulation in 
which fcrmionic degrees of freedom arc associated to both 
the scale factor and to the scalar field. 

By realizing the algebra of the fermionic variables and 
representing them as matrices, we get four equations for 
four components of the wave function. We focus our at- 
tention in two of them that are independent, and apply 
the WKB method in order to get two classical SUSY- 
cosmological equations. The associated equations of mo- 
tion for the scalar field are obtained by means of Hamil- 
ton's equations. 

In these supersymmetric Eisntein-Klein-Gordon equa- 
tions (SUSY-EKG), new contributions arise that behave 
like stiff matter, and some others in which the usual 
scalar field terms are modified by functions of the scale 
factor. 



For simplicity, we focused our attention in a flat Uni- 
verse and were able to find exact solutions of the equa- 
tions of motion. In one of them the scalar field potential 
is a negative constant, and then the radius of the Uni- 
verse is a periodic function. A second exact solution of 
interest corresponds to the case of an exponential (su- 
per)potential. This is a scaling solution corresponding to 
stiff fluid matter as revealed by the power law behaviour 
of the scale factor. 

We performed an analysis of the dynamical system 
structure of the SUSY-EKG equations in order to find all 
relevant physical solutions. Not surprisingly, we found, 
basically, the same solutions that appear in the stan- 
dard classical case. However, the stability and existence 
properties of the solutions were strongly modified by the 
supersymmetric corrections. 

In general, we can say that all solutions show that the 
scale factor and the scalar field degrees of freedom behave 
like in the case of stiff matter domination, because the 
inflationary solution is absent. This a consequence of the 
correction terms that appear in the equations of motion, 
in which the scalar field functions are mediated by stiff- 
matter terms of the scale factor. 

All conclusions are consequences only of the supersym- 
metric nature of the equations of motion, and the solu- 
tions are expected to be relevant in the early stages of the 
Universe. We cannot foretell the consequences that may 
arise in the case of a less exact supersymmetry, or even 
a broken one. But these are possibilities that may mod- 
ify again the phase structure of the solutions and allow 
the existence of inflationary solutions. This is research 
beyond the purposes of the present work that we expect 
to report elsewhere. 
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